The effects of sub-shells in highly magnetized relativistic flows 



Jonathan Granot 



ABSTRACT 

Astrophysical sources of relativistic jets or outflows, such as gamma-ray bursts 
(GRBs), active galactic nuclei (AGN) or micro-quasars, often show strong time 
variability. Despite such impulsive behavior, most models of these sources assume 
a steady state for simplicity Here I consider a time- dependent outflow that is 
initially highly magnetized and divided into many well-separated sub-shells, as 
it experiences impulsive magnetic acceleration and interacts with the external 
medium. In AGN the deceleration by the external medium is usually unimportant 
and most of the initial magnetic energy is naturally converted into kinetic energy, 
leading to efficient dissipation in internal shocks as the sub-shells collide. Such 
efficient low-magnetization internal shocks can also naturally occur in GRBs, 
where the deceleration by the external medium can be important. A strong low- 
magnetization reverse shock can develop, and the initial division into sub-shells 
allows it to be relativistic and its emission to peak on the timescale of the prompt 
GRB duration (which is not possible for a single shell). Sub-shells also enable the 
outflow to reach much higher Lorentz factors that help satisfy existing constraints 
on GRBs from intrinsic pair opacity and from the afterglow onset time. 

Subject headings: gamma-rays burst: general — magnetohydrodynamics (MHD) 
— shock waves — ISM: jets and outflows 



1. Introduction 

The composition and acceleration mechanism of the relativistic outflows that power 



$amma-ray bursts (GRBs) are important open questions in this field (for a review see 
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20051 ). In particular, their degree of magnetization and the role of magnetic fields in the 
acceleration or collimation of GRB outflows is of great interest. In recent years, models in 
which the outflow is highly magnetized close to the source, and possibly also at very large 
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m & Spruit 
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Vlahakis & Koniel 


2003; 
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tov & Blandford 


2003; 
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2006; 
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200ft 


Tchekhovskov et al. 


2010; 


Lvubarskv 


20101 


)) and may provide a viable alternative to the traditional fireball model. 



Some other relativistic outflow sources are even more likely strongly magnetized near the 
central source. Pulsar winds are almost certainly Poynting flux dominated near the source, 
and the same very likely also holds for active galactic nuclei (AGN) and tidal disruption 
events (TDEs) of a star by a super-massive black hole. In AGN and TDEs, which are often 
highly variable (i.e. impulsive) suggesting sub-shells in the outflow, since the central accreting 
black hole is super-massive then even close to it the Thompson optical depth tt may not 
be high enough for thermal acceleration by radi ation pressure - the main competition to 
magnetic acceleration - to work efficiently (e.g., iGhiseUinil 1201 ll ). Observations of relevant 
sources, such AGN, GRBs or pulsar wind nebulae suggest that the outflow magnetization is 
rather low at large distances from the source. An important outstanding question concerning 
outflows that start off highly magnetized near the source is how they convert most of their 
initial electromagnetic energy into the energy in the bulk and random motions particles, 
where the latter also produces the radiation we observe from these sources. This is known 
as the a problem, namely how to transform from a ^> 1 near the source to o <C 1 very far 
from the source, where the magnetization parameter a is the Poynting-to-matter energy flux 
ratio. Different approaches to this problem have been considered so far. 

Outflows that are initially Poynting flux dominated are usually treated (for simplicity) 
under ideal MHD, axi-symmetry and steady-state. Under these conditions, however, it is 
difficult to achieve sufficiently low magnetization io < 1 or even o <C 1) at large distances 



from the source that would a llow efficient dissipation in internal shocks (IKomissarov et al. 



2009; 


Lvubarskv 


2009. 


2010a) 



(a ^> 1) also at large distances from the source and the observed emission is powered by 



magn etic reconnection rather than by internal shocks (ILyutikov fc Blandfordll2003l ; iLyutikov 



20061 ). Alternatively, the non-axi-symmetric kink instability could randomize the direction of 



the magnetic field, making it behave more like a fluid and enhancing ma gnetic reconnection. 



which both increase the acceleration and help lower the m agnetization (IHeinz fc Begelman 



2000 ; iDrenkhahn fc Spruitl 120021; iGiannios fc Spruitl 120061 ). Another option that may be 



relevant for AGN and GRBs (ILyubarskyi l2010bl ). is that if the Poynting flux dominated 



outflow has alternating fields (e.g. a striped wind) then the Kruskal-Schwarzschild instability 
(i.e. the magnetic version of the Rayleigh- Taylor instability) of the current sheets could lead 
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to significant magnetic reconnection, which in turn increases the initial acceleration resulting 
in a positive feedback and self-sustained acceleration that leads to a low magnetization. 

Here I focus on the effects of strong time dependenc e - impulsive outflows that are ini- 



tially highly magnetized, with <jq 3> 1, under ideal MHD. iGranot. Komissarov fc Spitkovsky 



( 1201 ll . hereafter paper I) have recently found a new impulsive magnetic acceleration mech- 



anis m for relativ i stic o utflows, which is qualitatively different from its Newtonian ana- 



log ( IContopoulosI Il995l ). and can lead to kinetic energy dominance and low a that allow 
for efficient dissipation in internal shocks. Paper I focused mainly on the acceleration of an 
initially highly magnetized shell of plasma into vacuum. The initial magnetic energy can be 
almost fully converted into kinetic form as the shell expands radially under its own magnetic 
pressure. Initially, while it is still highly magnetized, this expansion leads to impulsive ac- 
celeration and the resulting increase in the Lorentz contraction almost exactly cancels the 
increase in the shell's width in its own rest frame, leading to a constant width in the lab (i.e. 
central source) frame. Once it becomes kinetically dominated it starts spreading radially 
significantly also in the lab frame, and its magnetization quickly drops to a <C 1. 

Paper I only briefly discussed the effects of the interaction with the external medium. 
The interaction with an unmagnetized external medium whose density varies as a power- 
law wit h the distan ce from the central source is analyzed in detail in an accompanying 



paper ( jGranotll20! ll hereafter paper II). The present work generalizes this self-consistent 
treatment of the combined impulsive magnetic acceleration and deceleration by the external 
medium, by examining the effects of an outflow that is initially divided into many well 
separated sub-shells, instead of a single shell. Such multiple sub-shells are both naturally 
expected in highly variable sources, and are also required in order to produce internal shocks 
(where enabling efficient internal shocks is one of the main motivations for this type of 
model). 

The basic test case examined in paper I was a highly magnetized shell initially at rest, 
whose back leans against a conducting wall and whose front edge faces vacuum, while in paper 
II the vacuum is replaced by an unmagnetized external medium. As discussed in paper I, 
while such a setup might be directly applicable for, e.g, a giant flare from a soft gamma 
repeater, for most relevant astrophysical sources (such as GRBs, AGN or micro-quasars) we 
expect that initially a quasi-steady acceleration takes place and saturates at some radius (i.e. 
distance from the source), while the impulsive acceleration mechanism takes over at a larger 
radius. The dynamics most relevant for our purposes occur after the impulsive acceleration 
takes over, and are insensitive to whether part of the earlier acceleration occurred in the 
quasi-steady regime (while the jet was being collimated). However, during the impulsive 
acceleration phase the outflow still retains memory of its initial properties when it was 
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ejected from the central source, namely the time history of its initial magnetization o"o and 
ejected energy per unit time, which is reflected in the isotropic equivalent luminosity L (i.e. 
energy flux through a sphere of fixed radius if the outflow occupied all of the solid angle). 
During this phase the outflow opening angle remains roughly constant, and its dynamics are 
essentially spherical, and equivalent to those for the planar shown in paper I. 

The sub-shells are assumed to be initially well separated, uniform with sharp edges 
and separations comparable to or larger than their widths, and with a large contrast (e.g. 
in energy density) between the sub-shells and the "gaps" between them (as might be the 
result of large variations in the energy output rate from the central source into the outflow). 
A modest contrast and/or very smooth edges for the sub-shells might cause the sub-shells 
to collide and merge earlier on, while they are still highly magnetized, making the outflow 
subsequently behave closer to the steady regime than to the impulsive regime. Quantifying 
these effects, however, is beyond the scope of this paper and is left for a separate work. 

The basic physical setup explored in this work is described in § [2j while § [3] studies 
in detail the main test case, namely the dynamics of iV 3> 1 identical sub-shells. I find 
several different cases for the dynamics, which naturally divide into three groups: (i) a low 
magnetization "thin shell" (cases 1 and 2*) where a strong mildly relativistic reverse shock 
develops on a timescale larger than the duration of the prompt GRB emission (Td ec ^> ^grb), 
(ii) a low or mild magnetization "thick shell" (cases 2 and 3A), where a relativistic reverse 
shock develops whose emission peaks on a timescale comparable to that of the prompt GRB 
emission {T^ ec ~ Tqrb), (in) a high magnetization "thick shell" (case 3B), where the reverse 
shock and its associated emission are strongly suppressed (with T& ec ~ Tqrb where T& ec in 
this case can be identified through the afterglow emission, which peaks on this timescale). 
The first two groups, (i) and (ii), naturally produce internal shocks at mild magnetization 
and high Lorentz factors, which can help accommodate GRB observations. The effects of 
varying the initial magnetization between different sub-shells are explored in § HJ and I show 
that this is not expected to strongly affect the main results, even for strong initial variations 
(and the same holds for order unity variations in the other model parameters). The new 
results found in this work are summarized in § [5] and their implications are discussed in § [6j 



2. Sub-shells versus a single shell 

Here the results of paper II for a single spherical shell of initial (lab-frame) width 
A ~ Ro, energy E, luminosity L « Ec/A and initial magnetization cr = Bq / (47r pod 2 ) ^> 
1 (with a magnetic field normal to the radial direction, where Bq and po are its initial 
magnetic field and rest-mass density, respectively), are generalized to the case where the 
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same initially uniform shell is divided into iV 3> 1 identical sub-shells of initial widths 
A 0i sh = A /iV and separations A gap > A 0iS h, energy E sh = E/N, with the same luminosity 
(L sh = E sh c/A 0jSh = Ec/A = L) and initial magnetization (o" 0iS h = °o) as the original 
single shell. The total initial width of all the sub-shells including the gaps between them is 
A tot « (A 0iSh + A gap )iV = A A where A = 1 + A gap /A , s h- 

The unperturbed external medium is taken to be cold, unmagnetized, and with a rest 
mass density that varies as a power-law (of index k) with the distance R from the central 
source, p\ = AR~ k . The external medium interacts directly only with the first (i.e. leading) 
sub-shell, which sweeps it up and drives a strong relativistic shock into it, where a contact 
discontinuity (CD) separates the shocked external medium and the material in the first 
sub-shell. The subsequent (or trailing) sub-shells do not directly interact with the external 
medium. Instead, they interact directly only with their neighboring sub-shells (i.e., each sub- 
shell interacts only with the sub- shell/s just in front of it and/or just behind it), and initially 
they propagate in the relatively evacuated region left behind by the preceding sub-shell. 

Since the sub-shells are assumed to be initially well separated, each sub-shell starts ac- 
celerating independently. The first (or leading) sub-shell interacts with the external medium, 
and until the second sub-shell collides with it from behind its dynamics essentially follow 
those of a single (albeit, relatively narrow) shell, which were studied in detail in paper II. 
Following the results of paper II (and adding 'sh' to the subscripts of quantities that refer 
to a single sub-shell rather than to the whole outflow), it has the following critical radii: 

Ro* ~ ^ , i2c>h ~ § , ^,sh ~ R U N-^ 10 -^ , R cr , sh ~ i? cr 7V- 2 /(4-*) . (1) 

Here i?o,sh (Ro) is the initial radius of the sub-shell (or the original single shell), essentially 
equal to its initial width, A 0iS h (Ao), where its typical Lorentz factor and magnetization 
are (r sh ) ~ a^ h ((V) ~ o-J /3 ) and (cr sh ) ~ <rj^ ((a) ~ (X 2/3 ), respectively, and in regimes 
I or II from paper II at R > i?o,sh (R > Ro) they start to evolve as (r s h) ~ 0"o/(°"sh) ~ 
(o- ,sh-R/-Ro,sh) 1/3 ((r) ~ cr /(a) ~ (o- ^/-Ro) 1/3 ); -R c ,sh (Rc) is the coasting radius where 
if the external density is low enough, corresponding to regime I from paper II, the sub- 
shell (or the original single shell) becomes kinetically dominated and starts to coast at 
(r s h) ~ o"o )S h ((r) ~ o"o) while its magnetization rapidly drops with radius, (cr s h)(-R > 
Rc,sh) ~ R c ,sh/R > Rc) ~ Rc/R)', R u ,sh (R u ) is the radius where, in regime II from 

paper II, the typical Lorentz factor of the sub-shell (or the original single shell) becomes 
similar to that just behind the CD, stops growing as R 1 ^ 3 , and instead starts evolving as 
(r sh ) ~ a /(a sh ) ~ r criSh (i?/ J R cr , sh )( fe - 2 )/ 4 ((r) ~ a /(a) ~ T^R/R^-^) up to the 
radius -R C r,sh (Rcr)', -Rcr.sh ~ -Ro,shrcr,sh (-^cr ~ Ro^cr) * s the deceleration radius in regime II 
(from paper II) where most of the energy originally in the sub-shell (or the original single 
shell) is transfered to the shocked external medium. 
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l°"0,sh 



The critical Lorentz factors, r cr s h or r cr , which signify the borderline between regimes I 
< T crsh or a < IV - A TT {r r( 12 " 3fe )/ 2 ~r ^ ^ t^-»>/^ 



and II (r criSh < a 0iS h < r^ z sh OKJ/z or T cr < a < 0,ej/z ), satisfy 



rcr ' sh ^ N (2-k)/(8-2k) 



R 



cr,sh 



R r 



(fc-2)/4 



(2) 



This essentially reflects the fact that since the sub-shell and the corresponding (or orig- 
inal) single shell have the same luminosity, they share the same line in the T — R plane 
corresponding to a pressure balance at the CD between the bulk of the shell (or sub-shell) 
and the shocked external medium (for details see paper II), Tcd ~ [L/A(?) 1 ^R t - k ~ 2 ^ A ~ 



J- cr (R/Rcr) 



(fc-2)/4 



J- cr,sh 

(R/R 

cr,sh 



,(*-2)/4 



When comparing to observational constraints, 



the numerical value of r cr is very useful, 



(3 - k)E 
4nAc 2 A 3 - k 



(3-fc)(l + z) 3 - k A 3 - k E 
A * Ac5 ~ kT GRB 



395 ( 3 / 8 A 3 / 8 EH 8 n 1/8 T- 3/8 (k = 0) 



88 C 1/4 A 1 /4 j e; 5 1 3 /4 A; 1/4 T3- 1/4 (k = 2) 



(3) 



(paper II) where ( = (l + z)/3, z is the source redshift, Tgrb = (1 + z)A tot /c = 30T 30 s is the 
observed duration of the GRB and E = 10 53 i?53 erg is the total (isotropic equivalent) energy 
of the ejecta. Numerical values are provided for the physically interesting cases of k = (a 
uniform medium of number density n = A/m p = n cm -3 ), and k = 2 (corresponding to the 
stellar wind of a massive star progenitor, with A = 5 x 10 n A* gr cm -1 ). The subsequent 
sections largely follow paper II, build upon it and use consistent notations. 



3. Identical sub-shells: dynamical regimes for different values of k and <Jo iS h 

3.1. Moderate (or no) external density stratification: k < 2 

Case 1: If cTo jS h = cro < Tcr (or -R CiS h < R cr /N; top panel of Fig. [TJ) then even the 
first (leading) sub-shell reaches its coasting radius -R c>s h without being significantly affected 
by the external medium. Since the subsequent (trailing) sub-shells start propagating in 
the evacuated region behind the first shell, they are initially not affected by the exter- 
nal medium. Thus, all of the sub-shells accelerate largely independently until becom- 
ing kinetically dominated at -R c ,sh ~ -Ro,shC"osh ~ Rc/N, where they start coasting (at 
(r s h) ~ °"o,sh) and spreading radially. Then, they soon collide with each other and form 
internal shocks at Ris ~ (A gap /A 0jSh )i? csh . After the internal shocks subside, the resulting 
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merged shell of width ~ A to t ~ A A ~ A (A gap /Ao,sh) starts spreading radially only at 
R s ~ A tot a 2 sh ~ (Ag a p/A 0i sh)-Rc ~ NRis where its mean magnetization is still similar to 
that near the internal shocks radius, (a)(R s ) ~ (a)(Ri S ) ~ A 0iS h/A gap < 1. At R > R s 
it starts evolving as A ~ (R/R s )A tot ~ (R/R C )A and (a) ~ (R s /R)(cr)(R s ) ~ R c /R, 
similar to the single wide shell that corresponds to the sub-shells considered here if they 
had no initial gaps between them. Thus, the merged shell similarly follows the unmag- 
netized "thin shell" case (regime I from paper II), and the reverse shock that forms be- 
comes mildly relativistic near the deceleration radius, i?dec ~ -Rr ~ (E / OqAc 2 ) 1 ^^^ ~ 
-Rc(rcr/o"o) ( ' 8_2fc ' ) ^ 3_/c ' ) (here Rr is the radius at which a rest mass E /T^c 2 of the external 
medium is swept up; see Eq. [4] of paper II), where it finishes crossing the merged shell 
and most of the energy is transfered to the shocked external medium. The magnetiza- 
tion at this radius is already low, (a)(R dPr ) ~ fcrn/T rr ) (8 ~ 2fc) / ( 3 ~ fc) ~ Rj R dec < I. At 
R > R^ec the flow quickly approaches the iBlandford &: McKed (119761 ) self-similar solution, 



where (r) ~ (E/Ac 2 ) 1/2 R^ k 3)/2 ~ a (R/R r Y k 3)/2 . The maximal Lorentz factor reached by 
the outflow in case 1 is T ~ o"o < r cr . 

Case 2: If T cr < a = a , sh < r cr , sh (i.e. 1 < or /r cr < ^(2-^/(8-2/=)) or Rc j N < ^ < 
Rcr,sh ~ R CT N~ 2 ^ 4 ~ kS> (middle panel of Fig. [T]), then similarly to case 1 above all of the 
sub-shells accelerate independently until becoming kinetically dominated at -R CjS h, where 
they start to coast (at (T s h) ~ cr 0)S h) and spread radially, and soon collide with each other 
at Ris ~ (A gap /Ao iS h)-Rc,sh- The main difference is that in this case the resulting merged 
shell is decelerated significantly by the external medium before it starts to spread radiallyEl, 
while its typical magnetization is modest, ~ (a) (Ris) ~ A 0sh /A gap < 1, and therefore 
a strong, highly relativistic reverse shock develops. Hence, case 2 effectively reverts to the 
unmagnetized (or at most mildly magnetized) "thick shell" case, where a bright reverse shock 
emission on a timescale comparable to that of the prompt GRB emission (Tjg,, ~ ^grb) is 
expected. Note that both the observed deceleration time, T dec ~ R dec /c{T) 2 (R dec ) (i.e. the 
timescale on which both the reverse shock and afterglow emission peak) and Tgrb scale 
linearly with A in this case so that their ratio (or the relation T dec ~ Tgrb) is independent 
of A. The effective luminosity of the merged shell is somewhat lower than that of the 
original sub-shells or the corresponding single wide shell, L merged w Ec/A tot = L/A, and 
therefore a strong relativistic reverse shock develops at the radius R2 where the Lorentz 
factor of the CD at larger radii, T C b(R2 < R < Rdcc) ~ (L mcrgcd /Ac 3 ) 1/4 R^ k ' 2 ^ A , becomes 
comparable to the coasting Lorentz factor, Tcd(R < R2) ~ (r)(-R c ,sh < R < Rdec) ~ o"o,sh- 
This implies R2 ~ A _1 ^ 2_fe ^i?i wherd^l R\ is the radius at which a = 1 just behind the CD 



^his follows from the relations R s /R 2 > R s /Rdcc - K^-^'^-^ R c / R a - A( 3 - fe )/( 4 - fc ) (cr /r cr ) 2 > 1. 
2 The reason why R 2 is so close to Ri, up to a factor of ~ (i m orgcd/^) 1 ^ 2_fe ' ~ ^- 1 /( 2 ~ k ) i [ s 



-8 - 



(for the corresponding single wide shell) and it is given by R\/R c ~ {T cr / oq)^ 2 ^^ 2 ^^ (see 
Eqs. [27] and [34] in paper II). The relativistic reverse shock finishes crossing the shell at 
the deceleration radius, i?dec ~ A 1 ^ 4- ^^!., where most of the energy is transfered to the 
shocked external medium. This is the reason behind the sharp drop in the energy weighted 
mean Lorentz factor of the flow, (r), near Rdec (see middle panel of Fig. [1]) since essentially 
on a single dynamical time it changes from being dominated by the coasting (unshocked) 
part of the merged shell with V ~ o"o jS h to being dominated by the shocked external medium 

with r(i? dec ) ~ r cr A-(3-fc)/(8-2*) < 

°"o,sh (where near the transition, at R ~ -Rdec, the 
shocked part of the merged shell, which has a Lorentz factor similar to that of the shocked 
external mediu m, also holds a good f r action of the total energy). At R > Rd ec the flow quickly 
approaches the iBlandford &: McKed (119761 ) self-similar solution. The maximal Lorentz factor 
reached by the outflow in case 2 is (r) ~ o"o, and can approach r cr s h ~ r^A^ 2- ^^ 8-2 ^ > T cr . 

Case 3: If cT , s h = > r cr , sh (i.e. a /r cr > jV^)/^)) or R c ^ > R cv , sh ~ R CV N~ 2 /^ 
then the first sub-shell starts to decelerate significantly (because of the PdV work that 
it performs across the CD on the shocked external medium) at R u , s h ~ R u N^ 4 ^ w ^ 3k ^ ~ 
#cr, s h(r cr , s h/a , s h) 4/(10 ~ 3fc) , while it is still highly magnetized ((a sh )(R U:Sh ) ~ cr , s h/^ci){Ru, s h) ^ 
( cr o,sh/r cr>s h) ( - 8 ~ 2A '' )//< ' 10_3A '' ) > 1) and before spreading radially appreciably, thus effectively fol- 
lowing regime II (from paper II). The second sub-shell, however, would effectively collide 
and merge with the first sub-shell only at a radius -R co i,i ~ (Ag ap /Ao i sh) 2 ' /( - 4 ~ fc ' > -Rcr,sh (see 
Appendix 12} , and until that radius the bulk of this sub-shells would still accelerate almost 
as if into vacuum, following (r sh ) ~ cr 0iS h/( cr sh) ~ (c r o,sh-R/-Ro,sh) 1,/3 - The same holds for 
subsequent collisions, where the radius of the n th collision is -R co i,n ~ {nK/N) 2 ^ A ~ k ' > R cr ~ 
{n/N) 2 ^ A ~ k ^R& ec (see Appendix [A]), as long as -R co i,n < -R c ,sh or equivalenthyB n/iV < (cTo/r 3 B) 4 ~ 
(which is always satisfied in case 3B below, but not in case 3A), where is defined be- 



that in both cases it is basically the radius where the Lorentz factor of the contact discontinuity (CD), 
Tcd ~ (L/Ac 3 ) 1 / 4 i?( fc ~ 2 - ) / 4 , which is determined by pressure balance at the CD between the shocked exter- 
nal medium and the bulk of the shell (where the latter is provided by magnetic pressure for R\ in regime II 
and by thermal pressure behind the mildly relativistic reverse shock for -R2), is ~ 00 = <T 0iS h- 

3 The front part of the sub-shell would start interacting with the tail of the preceding sub-shell at smaller 
radii, but the interaction would start to significantly affect the bulk of the second sub-shell only at R ~ i? C oi,i- 
By i? C oi,i the first sub-shell has already expanded significantly in the lab frame, by a factor of ~ A, and 
therefore its electromagnetic energy is reduced by the same factor (since it scales as -Eem °c B 2 A oc 1/A), 
however its magnetization is still high so that most of its energy is still in magnetic form and this represents 
mainly the loss of energy to the shocked external medium due to the work it performs on it at the CD. 

4 If R co \, n > Rc,sh or n/N > (cr /r 3 B) 4_ ' c , then the sub-shell would spread radially and collide at Ris ~ 
Ai? c . s h, and would eventually be decelerated by the relativistic reverse shock at R ~ R co i,n (which might in 
that sense still be considered as an effective "collision" radius). 
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low. In such highly magnetized collisions the (n + l) th sub-shell that catches up from be- 
hind in the n th collision accelerates up to the time of that collision to (r sh)ri+1 ) (R co i,a) ~ 
(d , sh R cohn /R ,s h ) 1/3 ~ <ToARcoi,jR c , aii ) 1/3 ~ r 3 B(ao, s h/r 3 B) 1/3 (^/iV) 2/(12 - 3fc) ,'and reaches a 
magnetization (a ah>n+1 ) (R col>n ) ~ (i? co i,n/^c,sh)- 1/3 ~ Ksh/r 3B ) 2/3 (n/A0- 2/(12 - 3fe) . 

Case 3A: for R CT , sh < -R CjSh < -R dC c ~ -R cr A 1/(4 - fe) (or r crjSh < cx ,sh < r 3A where T 3A ~ 
(Rdec/ Ro,shY^ 2 ~ T^N^-^A 1 ^ 8 ^ 2 ^ is the maximal attainable Lorentz factor in the shell with 
mild magnetization, <j < 1; bottom panel of Fig. [T]) the trailing sub-shells eventually (for 
n/N > (cr /r 3B ) 4_ ' c ) become kinetically dominated at -R CjS h- They then start coasting and 
spreading radially, thus quickly colliding and merging at Ri$ ~ (A gap /Ao iS h)-Rc,sh where 
(a) ~ A 0i sh/Ag ap < 1, allowing for efficient energy dissipation in the resulting modest 
magnetization internal shocks. The resulting merged modestly magnetized shell is then 
gradually decelerated by a relativistic reverse shock, largely following the unmagnetized 
"thick shell" case at R > Ris (so that i?dec ~ -RcrA 1 ^ 4- ^, similar to case 2). Since -R CiS h < 
Rdec corresponds to cr 0sh < T 3 a, this implies that in case 3A the Lorentz factor of the 
resulting internal shocks can be as high as (r) ~ <7o < T 3 a ~ r cr A^ 1 / 2 A 1 ^ 8_2A ''\ or up to a 
factor of ~ A^ 1 / 2 A 1 /( 8 - 2fc ) > 1 larger than T cr . 

Case 3B: for i? c , sh > R dcc ~ i? cr A 2 /( 4 ~ fc ) (or o- , sh > T 3B ~ T^N^A 1 /^) even the last 
sub-shells are still highly magnetized and do not spread radially appreciably by the time 
they collide and merge with the preceding sub-shell (as the latter is decelerated by the ex- 
ternal medium or the preceding sub-shell and spreads radially just before colliding with the 
subsequent sub-shell). Thus, subsequent collisions with later ejected sub-shells proceed at 
the back end of the growing, highly magnetized, quasi-uniform, merged shell behind the 
CD, whose Lorentz factor evolves (on average) as r merge d ~ (-^merged/ ' Ac , ) l / 4 R <yk ~ 2 ^ i , until 
the deceleration radiu^l, R& ec ~ R cr K 2 ^ A ~ k \ The last collision occurs at R co \,n ~ -Rdec and 
the last shell reaches the maximal Lorentz factor of (r s h,jv) ~ r 3 B(tT 0jS h/r 3 B) 1//3 , which is 
larger than r 3B by a factor of (cro^/T^) 1 / 3 > 1. Since -R co i,n ~ Rdedn/N) 2 ^ 4 '^, most 
of the collisions occur rather close to -Rdec- Hence, case 3B largely follows the highly- 
magnetized "thick shell" case of regime II, with L — > L merge d ~ L/A 2 , where (T)(R < 
-Rdec) ~ (°~o,shR/ Ro,shY^ 3 is dominated by the almost freely expanding and accelerating 
shells at the back of the flow. Near -Rdec; essentially within a single dynamical time, all 



5 In case 3B i?dec is slightly larger than in cases 2 or 3A because the sub-shells remain highly magnetized 
near i?dec oc Lmergcd^dcc) -1 ^ 4 "^ , so that the energy of each sub-shell decreases by a factor of ~ A while 
its width grows by a similar factor, resulting in L roer ged ~ L/A 2 and i?dec ~ i? cr A 2 /( 4-fc ). In case 2 (or 
3A), all (or late) sub-shells become kinetically dominated and their energy remains constant until they are 
decelerated by the reverse shock, while their width increases by a factor of A well before -Rdcc, so that 

imcrgcd ~ L/A and i? dcc ~ i? CI A 1 /(4-fc). 
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of the remaining freely expanding shells collide and merge, so that (r) moves from be- 
ing dominated by those shells with a typical Lorentz factor (r)(i? dcc _) ~ F^ao/T^) 1 / 3 
and magnetization (a)(Rdec~ ) ~ cr o/(r)(-Rdcc,-) ~ (cr /T- iB ) 2 / 3 to being dominated by 
the unmagnetized shocked external medium with (r)(i?dec,+) ~ T CT A^^ k ^^ k \ where 
near i?dec a significant fraction of the total energy also resides in the merged highly mag- 
netized outflow shell, which has a typical Lorentz factor (r)(i?dec,+) and magnetization 
(cr)(R dec ,+) ~ a /A(T)(R dec>+ ) ~ A-^-kWl^ » i. 



3.2. Stronger external density stratification: 2 < k < 3 

Case 1 now corresponds to cr 0iSh = cr < r criSh (i.e. cr /T cr ~ (R c /R cr ) 1/2 < AT( 2 - fc )/( 8 - 2fe )) or 
i? CjS h < -R cr ,sh (i-e. -R CjS h/-Rcr < iV -2 ^ 4 "^ ) , but otherwise the behavior of the system in this 
regime is very similar to that in case 1 for k < 2 (see top panel of Fig. [2]). 

Case 2*: what used to be case 2 for k < 2 now corresponds to r^sh < er = cro.sh < r cr 
(i.e. iV( 2 - fe )/( 8 - 2fc ) < a /r cr < 1) or i? cr , sh < i? c , sh < i? , sh r 2 r (i.e. N~ 2 /^ < R c ', sh /R CT < 
N^ 1 ), and is called case 2* since its properties are somewhat different (see middle panel 
of Fig. [2]). The first sub-shell is in regime II, i.e. it starts to be significantly affected 
by the external medium (accelerates more slowly, as (r s h) oc _R( fc-2 )/ 4 instead of _R 1//3 ) at 
R u ,sh ~ RuN'^t 10-3 ^ . The first few subsequent sub-shells would still collide with the back 
of the growing highly magnetized merged shell behind the CD. However, once the radius 
of this merged shell (or of the CD) exceeds -R c ,sh the sub-shells first become kinetically 
dominated at -R c>s h, start coasting at (r sh ) ~ Oo )S h, expand radially, and collide with each 
other (at Ris ~ (A gap /A sh )_R C)Sh ) with a moderate magnetization ((c)(i?is) ~ A 0) sh/A g ap ^ 
1) before effectively colliding with (or being decelerated by) the back end of the merged 
shell. Therefore, once they start being decelerated by the merged shell behind the CD 
(which has F ~ Tcd) it occurs in the form of a highly relativistic reverse shock, as long as 
Tcd *C o"o,sh- However, this shock becomes Newtonian and weak at R2 ~ A _1 ^ 2_fc ^i?! (where 
Tcd ~ ^sh)) an d from that point on the merged mildly magnetized shell essentially coasts 
at T ~ cr ,sh = cr . The resulting merged shell of width ~ A tot w A A ~ Ao(A gap /Ao, s h) 
starts spreading radially only at R s ~ A tot er 2 sh ~ (A gap /Ao iS h)-Rc ~ NRjs where its mean 
magnetization is still similar to that near the internal shock radius, (a)(R s ) ~ (a)(i?is) ~ 
A 0jSh /A gap < 1, and evolves similarly to case 1. At R > R s it starts evolving as A ~ 
(R/R s )A tot ~ (R/R C )A and (a) ~ (R g /R)(a)(R g ) ~ R c j R, similar to the original single 
wide shell, thus following the unmagnetized "thin shell" case. One possible difference is that 
the early short-lived phase of a strong relativistic reverse shock might result in an observable 
early and short-lived spike in the reverse shock emission, on a timescale of ~ (i?2/-Rr)7dec ~ 
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(r cr /o"o)^ 8_2fc ^' ( - 3_ ' c ^ 2_A: ' ) 'Tdec that is precedes the main reverse shock emission peak (which 
peaks on a larger timescale of T dcc ). For k < 3 the shell is eventually decelerated at i?dec ~ Rr 
by the reverse shock, which becomes mild ly relativistic by that radius, and then the flow 
approaches the iBlandford &: McKed (119761 ) self-similar solution. 



Case 3 now corresponds to o"o jS h = cro > Tcr (i-e. cr /r cr ~ (R c / R cr ) l l 2 > 1) or -R CjS h > 
Ro,shT 2 cr ~ Rcv/N. Case 3A corresponds to Ro, sh T 2 cv < R c , sb < R dec (i-e. N' 1 < R c , sh /R CI < 
A 1 /(4-fc)) or r cr < a < r 3A ~ r cr A^ 1 / 2 A 1 /( 8 - 2fe ) (see bottom panel of Fig. E), while Case 
3B: corresponds to -R CiS h > -Rdec (i-e. -R CjS h/-R C r > A 2 /( 4 ~ fc )) or cr > r 3B . Other than that, 
cases 3A and 3B behave very similarly to k < 2 (described above). 



3.3. A wind-like external density profile: k = 2 

For k = 2 we have r cr sh = r cr and -R C r,sh = Ro,shTl r , so that there is no case 2 or 2*. 
Case 1 corresponds to cr ,sh = 0"o < r cr = r crsh (or -R CjS h < -R C r,sh = -Roarer)- Case 3 A 
corresponds to r criSh = r cr < cr ,sh < r 3A ) (or R CTtSh = Rq^T^ < R C)Sh < R dcc )- Case 3B 
corresponds to cr ,sh > r 3B (or R C)Sh > R dec ). 



4. Varying the initial magnetization <Jo iS h between different sub-shells 



It is reasonable to expect that the magnetization of different sub-shells might differ, at 
least by factors of order unity. Here we consider the effects of such a variation, while keeping 



the other sub-shell parameters fixed (namely A 



0,sh 



R 



O.sh- 



Agap; -^sh 



L). 



First, let us examine whether sub-shells might collide during the acceleration stage. 
Consider two sub-shells ejected with a time difference t gap ~ A gap /c, the first with o"o jS h = 
cr 0j i and the second with cr 0iS h = co,2 > Co,i- Each sub-shell initially accelerates as Tj ~ 
(cro,i-R/-Ro,sh) 1//3 and its width remains almost constant in the lab frame up to its coasting 
radius, R Ct i ~ i?o,shO"oi; so ^ na ^ both sub-shells accelerate at R < R Ci i < R c ^- Hence, 
when making the simplifying assumption of uniform sub-shells, the separation between them 
evolves as 



l(R < R 



c,l, 



Agap 



R 



A 



3A . sh 



dR 
IT 



1 



r 2 

1 2 



gap 



2a, 



2/3 
0,1 



O"0,2 



2/3' 



R 



R 



0,sh 



1/3 



(4) 



The two sub-shells would effectively collide when l(R) = 0, corresponding to a collision 
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radius R co \ 3> Ro, s h, which for R co \ < _R c ,i is given by 



3 



2/3 



3 




l - 




(5) 



Thus, for A gap > 1.5[1 - (o ,i/o , 2 ) 2 / 3 ] A 0j sh such a collision would occur at -R co i > R c \, after 
the first (and slower) shell reaches its coasting radius and becomes kinetically dominated. 
This condition always holds for A gap > 1.5A , s h, regardless of the ratio of the initial sub-sell 
magnetizations, 0-0,2/00,1 > 1, and reasonably low values of this ratio relax the condition on 
A ga p/A 0iS h (e.g. it becomes A gap /A , sh > 1 for o , 2 /o ,i = 3 3/2 « 5.2). 

For A gap /A 0i sh < 1 one cannot neglect the radial spreading of the first sub-shell before its 
coasting radius, i2 Cj i, since it would result in a collision at Ris ~ (A gap /A , s h) 3 -Rc,i < R c ,i- 
For A gap /Ao iS h > 1 the effect of this radial spreading on the decrease in the separation 
between the two sub-shells becomes comparable to or larger than that of the difference in 
their typical Lorentz factors at R > R c ,i- Thus, when accounting for both of these effects, for 
A gap /A 0i sh ^ 1-1.5 the two sub-shells effectively collide near i? IS ~ (A gap /A , s h)-Rc,i > R c ,i- 
For a large contrast in the initial magnetizations, 00,2/00,1 > (A gap /A 0iS h) 1//2 , the second sub- 
shell would still be highly magnetized at the collision radius, (o sh , 2 )(-Ris) ~ (-Rc^-Ris) 1 ^ 3 ~ 
(A gap /Ao,sh) _1 ^ 3 (oo,2/oo,i) 2 ^ 3 > 1, while for a mild initial magnetization contrast, 00,2/00,1 < 
(A gap /A ,sh) 1//2 , it would be mildly magnetized, (o sh , 2 )(-Ris) < 1- The first sub-shell would 
always be mildly magnetized in this regime, (o s h,i)(i?is) ~ A , s h/A gap < 1. This would allow 
for reasonably efficient dissipation in the resulting internal shocks. 

If the trailing sub-shell has a lower initial magnetization than the leading one (o ,2 < 
o ,i) then the (effective) separation between the sub-shells initially increases during the 
acceleration stage. Once the trailing sub-shell becomes kinetically dominated and starts 
coasting at i? c , 2 < i? c ,i it starts spreading radially, but even if its head moves at very close to 
the speed of light it would effectively collide with the bulk of the first sub-shell only after the 
latter starts spreading appreciably. Such a spreading of the first sub-shell can occur either if 
it reaches its coasting radius (R c ,i) and becomes kinetically dominated (in which case both 
sub-shells would have a mild or low magnetization when they collide), or alternatively if 
it is still highly magnetized but decelerates as it transfers a good part of its energy to the 
sub-shell in front of it (or the shocked external medium across the CD) through PdV work. 

Altogether, a reasonable spread in the initial magnetization of the sub-shells, of <5o"o, s h ~ 
°o,sh) would not have a very large effect on the overall dynamics or the efficiency of the 
resulting internal shocks. Even a very high contrast, of 5o , s h o , s h, is expected to typically 
affect the overall efficiency of the internal shocks only by a factor of order unity (since at 
least one of the sub-shells in each collision is expected to have mild or low magnetization). 



13 



Moreover, it is not even clear whether the overall efficiency would actually be decreased or 
increased, since while some of the collisions might occur where one of the colliding sub-shells 
has a higher magnetization (compared to the case Sa 0:S h <C Co,sh)j the relative Lorentz factor 
of the colliding shells can be higher, thus increasing the the efficiency of the dissipation in 
the low magnetization colliding sub-shell. Similarly, order unity variations in A 0iS h, A gap , or 
L s h = i£shc/Ao,sh are n °t expected to have a very large effect on the overall dynamics or on 
the efficiency of the resulting internal shocks or reverse shock. 

As for the reverse shock, if some of the sub-shells remain highly magnetized within the 
merged shell that accumulates behind the CD, then this could suppress the reverse shock as 
it passes these sub-shells, thus effectively causing its emission to turn on and off as it passes 
regions (corresponding to different original sub-shells) of low and high magnetization. Such a 
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5. Summary 



To summarize, when the outflow consists of a large number of well separated sub-shells 
there are three main regimes. One of them - the highly magnetized "thick shell" (case 3B; 
#csh > #dec ~ .RcrA 2 ^ 4 -^ or R c , sh /R CI > A 2/(4 ~ fc) ) is expected to result m a severe suppres- 
sion of the reverse shock and its associated emission. Nevertheless, the thick shell nature 
of this regime can be deduced from the fact that the onset time of the afterglow emission, 
T dec , is similar to the duration of the prompt GRB emission (T dec ~ Tqrb). Therefore, this 
leaves two regimes for which a reasonably bright reverse shock emission may occur. In the 
low magnetization "thick shell" case (i.e. cases 2 or 3A; N^ 1 < i? CjS h/-R C r < A 1 ^ 4- *^), 
which was described in this work, the mean magnetization of the merged shell at the 
time when its bulk is crossed by the relativistic reverse shock (near i? d ec ~ A 1 / < - 4_fc ^i? cr ) 
is (o")(i?dec) ~ (ff)(Ris) ~ A 0! sh/A gap < 1, i.e. it is expected to be less than unity 
but not by a large factor (0.1 — 0.3 < (cr) (-Rdec) < 1)- Thus, very low values of (cr) 
are not expected in this regime, whose main observational signature is that the afterglow 
emission and the reverse shock emission both peak on a timescale similar to the dura- 
tion of the prompt GRB emission (T dec ~ T GRB ~ (1 + z)A tot /c). The expectations of 
this regime (both in terms of T dec ~ T GRB and the value of (cr)(.R dec )) appear consistent 
with the bright pro mpt optical emission from GRB 990123, which had been attributed 



to the reverse shock (lAkerlof et al.l Il999l ; ISari &: Piranl Il999t iFan et al.l l2002t IZhang et al. 



20031 ; iNakar &: Piranl 120051 ) . In the low magnetization "thin shell" case (regime I; cases 1 
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or 2*; R C:Sh /R cr < iV -1 ), the deceleration time T dec that corresponds to the duration of 
the peak reverse shock and afterglow emission components is expected to be larger than 
the prompt GRB duration, T dec 3> T G rb, and the magnetization at the radius where 
most of the energy is dissipated in the reverse shock (near i? dcc ~ i?r) is expected to 
be (a)(R dec ) ~ (Vr c r) 2{4 ~ fc)/(3 - fc) ~ R c /R dcc ~ T GRB /(T dec A) < 1, i.e. a factor of 
~ Tdcc/TcRB ^ 1 smaller than in the low magnetization "thick shell" case. Thus, a clear 
prediction of this model is a positive linear correlation between T GRB /T dec and (cx)(i? dec ) 
when there is a bright reverse shock emission (i.e. when {<j)(R dcc ) < A 0i sh/Ag ap < 1). 

Moreover, this model also has predictions for the internal shocks, which could be tested 
against observations if the prompt GRB emission is indeed from such internal shocks. The 
internal shocks radius is given by 



^R c sh ~ 10 14 ( 2 ( cm , (6) 

A ,sh ' A 0iSh S Vl0 2 - 5 7 V 0.1s J 1 ; 



which satisfies the usual relation, i? IS ~ r 2 (i? IS ) cT var obs , where T var obs is the observed 
variability time in the prompt GRB lightcurve. For reference, 

( 9.3 x 10 16 aV 4 C-V 4 A-V4 n -V4 s i/4 T V4 cm (fc = Q) ? 

#cr ~ <^ (7) 

[ 5.3 x lO^a^C^A-^A^Eli 2 ^ 2 cm (jfe = 2) , 
(see Eqs. [23] and [45] of paper II for the definition of a and the derivation of _R cr ) and 

10-^-3^/^-1/^-1/^/4 m y {k = oh 

(8) 



R 



c.sh 



Rr 



This shows that the relevant regimes correspond to reasonable model parameters, and could 
potentially occur in different GRBs. The most uncertain parameter is the initial magnetiza- 
tion, cro,sh or (To, whose value can be estimated in the low magnetization regimes, where the 
outflow becomes kinetically dominated (with T ~ o"o; lower limits on T also serve as lower lim- 
its on o"o in this model, without requiring kinetic dominance). Current constraints from pair 
opacity in the prompt emission and from the onset of the afterglow suggest 10 2 < cr < 10 3 . 
In case 3B, which corresponds to R CtS h/R C r > A 2 /( 4_fc ) > 1, all of the sub-shells collide while 
they are still highly magnetized, which suppresses the internal shocks and their associated 
emission, making them unlikely to power the prompt GRB emission (which in this case might 
be alternatively powered by magnetic reconnection events in the highly magnetized outflow). 
In all other cases the internal shocks occur at mild or low magnetization, allowing them to 
be reasonably efficient and potentially power the prompt GRB emission. 
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Discussion 



The effects of sub-shells in an impulsive, initially highly magnetized relativistic outflow 
have been studied, and compared to the case of a single wide shell. It has been shown that 
if a single wide uniform outflow shell is divided into a large number of sub-shells with a 
reasonable initial contrast and spacing between them (A gap > A 0iS h) then it could reach a 
significantly higher Lorentz factor^ The leading sub-shells effectively clear the way for the 
subsequent sub-shells, allowing them to accelerate for a longer time without feeling the effects 
of the external medium (almost as if into vacuum), thus enabling them to reach a higher 
Lorentz factor. Moreover, internal shocks arise from collisions between different sub-shells, 
which naturally occur at relatively high Lorentz factors and at low magnetizations that are 
vital in order to have a reasonable energy dissipation efficiency in the internal shocks. 

A sufficiently high Lorentz factor is needed t o overcome the compactness problem and 



1995 



avoid excessive pair production within the sou r ce (IKrolik &: Pierlll99ll ; iFenimore et al.l 11993 
Woods fc Loebl 



Baring fc Harding! 119971 ; iLithwick fc Sarill200ll ). It has been recently 
argued ( Levinson 201ol ) that the interaction with the external medium might not enable an 



impulsive highly magnetized outflow in GRBs to accelerate up to sufficiently high Lorentz fac- 
tors, and in particular that its maximal achievable Lorentz factor is largely limited to T < r cr . 
This would pose a particularly severe problem for a stellar wind-like external medium (k = 2) 
for which typically r cr < 10 2 (see Eq. [3]). Recent high-energy observations by the Fermi 
Large Area Telescope (LAT) have set a lower limit of T > 10 3 for the e mitting region in a 



numb er of GRBs with a bright high-energy emission ( lAbdo et al.ll2009al Jbl; lAckermann et al. 



20101 ) using a simplified one-zone model. Howev er, a more detailed and realistic treatment 



shows that the limit is 



2008; Ackermann et al. 



ower by a factor of ~ 3 (IGranot. Cohen- Tanugi & do Couto e Silva 
20111). which wou ld correspond to T > 10 2 5 for the brightest Fermi 



LAT GRBs (see also lHascoet et al.ll201ll ). This might nevertheless still pose a problem for 
a single highly magnetized shell in a stellar- wind environment. The present work, however, 
shows that if it is divided into a large number of sub-shells, then its Lorentz factor T could 
exceed r cr by up to a factor of ~ r 3 A/T cr ~ ]V 1 / 2 A 1 ^ 8_2fc ' ^> 1. Moreover, most of the dissi- 
pation in internal shocks is expected to occur near the maximal Lorentz factor attained by the 
outflow. This would greatly help satisfy the lower limits on T from compactne ss arguments, 
or from the onset of the afterglow emission (usually around a few hundred; ISari &; Piran 



iLevinsonl (|2010h has argued that for A gap ~ A 0jS h the sub-shells would effectively collide an merge well 



before their coasting radius (by a factor of ~ cr^ h » 1), while they are still highly magnetized ((o"sh) ~ 

<7q^ 9 3> 1), and would thus have a very small effect on the outflow, and in particular would not help to 
increase its maximal Lorentz factor. However, this conclusion is wrong and arises due to an error in his 
Eq. (29), which results in an incorrect expression for the collision time or radius. 
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19991 : iNakar fc Piranll2005l : Molinari et al.ll2Q07l : IZou fc Piranll2010l : iGruber et al.lboilf ). also 



for a stellar wind environment. 

It has been found (in papers I and II) that for a single shell there are two main dynamical 
regimes: the low magnetization "thin shell" (regime I, where the shell becomes kinetically 
dominated, coasts and spreads radially, and a reverse shock develops that becomes mildly 
relativistic near the deceleration radius, Rd ec ~ Rr, and whose emission peeks at Td cc ^> 
Tgrb), and the high magnetization "thick shell" (regime II or III, where the shell remains 
highly magnetized without reaching a coasting stage or spreading radially, and the the reverse 
shock is suppressed along with its associated emission). A high magnetization "thin shell" 
or a low magnetization "thick shell" are not possible for a single initially highly magnetized 
shell. In this work, it has been shown that a low magnetization "thick shell" regime becomes 
possible (and occurs in cases 2 or 3A) if such an initially highly magnetized shell is divided 
into a large number of sub-shells with reasonable initial separations (A gap > A , s h)- This 
would allow a relativistic reverse shock with bright emission on a timescale comparable to 
that of the prompt GRB emission (Td cc ~ Tgrb)- Moreover, if there are large variations 
in the magnetization between different sub-shells, this might result in alternating regions 
of high and low magnetization through which the reverse shock passes, causing it and its 
associated emission to be alternately suppressed and revived, resulting in a variable reverse 
shock emission. Such a reverse shock emission may bear some temporal correlation to the 
prom pt emission from the internal shocks (somewhat analogous to the pure hydrodynamic 
case; INakar &: Piranl 12004 ) . as both are affected by the magnetization of the sub-shells, 
though some delay might be expected between corresponding features in the internal shocks 
and in the reverse shock. 



The author thanks A. Spitkovsky, Y. E. Lyubarsky, T. Piran, A. Levinson and S. S. 
Komissarov for useful comments on the manuscript. This research was supported by the 
ERC advanced research grant "GRBs" . 
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Fig. 1. — The evolution of the typical or energy- weighted mean Lorentz factor ((r); thick 
solid magenta line) and the Lorentz factor of the CD (Tcd; thick dashed- dotted green line) 
with radius, R, for N 3> 1 equal, initially highly magnetized (cr 0sh ^> 1) sub-shells, compared 
to a single uniform shell (thick dashed line; purple in the top panel corresponding to regime 
I, and dark green in the middle or bottom panels corresponding to regime II), with the same 
initial magnetization (oq = o"o,sh) and luminosity (or energy density), as well as the same 
total energy and net width (not counting the initial gaps between the sub-shells), for k < 2. 
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Fig. 2. — The same as Fig. [T]but for 2 < k < 3. 



-21 - 



A. Estimating the collision radius of sub-shells with the uniform region 

behind the CD 



I consider here sub-shells that are still highly magnetized and have not spread radially 
significantly before colliding with the previous sub-shell that is decelerated by its PdV work 
on the shocked external medium across the CD (or on the preceding sub-shell). Let us begin 
with the first such collision. The first sub-shell is significantly decelerated by the external 
medium at -R UjS h ~ R u N~ A ^ 10 ~ zk \ at which stage its separation from the subsequent shell 
is still close to its initial value, A gap . One way of estimating the collision radius is that the 
head of the second sub-shell travels faster than the tail of the second sub-shell so that the 
difference in their velocities in the lab frame is A/3 
separation changes with radius as 
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Since the sub-shell's Lorentz factor before it is significantly decelerated is r sh (i? < R u 

(O"0,sh-R/-Ro,sh) 1//3 , We have i? UiS h/rc D (i? UiS h) ~ Ao i sh(-Ru,sh/-Rc,sh) 1 ^ 3 ~ Ao, s h(-Rcr,sh/-Ru,sh 

or A gap r^ D (i? UiSh )//? UiSh w (A gap /Ao,sh)(-Rcr,sh/-Ru,sh) (4 " fc)/2 - The collision occurs when the 
separation reaches zero, and in the relevant regime the second term in the square brackets 
in Eq. (lAlj) can can neglected, giving a collision radius of 
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A comparable estimate of i? co i,i is obtained when asking at what radius the radial spreading 
of the first sub shell due to the dispersion in its Lorentz factor (5T sa ~ (r s h)) and deceleration 
because of the work it performs on the shocked external medium across the CD becomes 
large enough to bridge the initial gap, A gap , from the subsequent sub-shell. This is since 



taking l(R) 
the shell, 



in Eq. f ]Alj) is essentially equivalent to the requirement on the spreading of 
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The above considerations can be readily generalized in order to estimate when a point 
with an initial lag of A moving at T ^> T CD catches up with the CD. Simply replacing A gap 
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with A ps nA ,shA gives the radius of the n'th collision, 

2/(4-*;) 



/ A v2/(4-fc) fnAY /n\2/(4-*) 

i?col,n ~ ( ) ^cr,sh ~ (nA) 2 /* 4 "^.* ~ I — 1 i? cr ~ i? dec , 

(A4) 

keeping in mind that i?dcc ~ A 2 ^ 4_fe ^i? cr in this regime. 



